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Abstract 



A simplified version of the Wigner-transformed time-dependent Hartree- 
Fock-Bogoliubov equations, leading to a solvable model for finite systems 
of fcrmions with pairing correlations, is introduced. In this model, pairing 
correlations result in a coupling of the Vlasov-type equation for the nor- 
| mal phase-space density with that for the imaginary part of the anoma- 

, lous density. The effect of pairing correlations on the linear response of 

the system is studied for a finite one-dimensional system and an explicit 
expression for the correlated propagator is given. 

2 ; 1 Introduction 

C 

The collisionless Boltzmann equation, or Vlasov equation, is a useful tool for 
studying quantum systems in the semiclassical limit ([1], p. 553). 
' It is well known that the Vlasov kinetic equation can be derived from the 

• time-dependent Hartree-Fock equation 

ihdtp = [h, p] (1) 

by taking the Wigner transform of the density matrix p, however, a straightfor- 
ward application of this formalism does not give immediately the Vlasov equa- 
tion: the additional approximation of neglecting the momentum dependence of 
the Wigner-transformed Hartree-Fock potential is required. 

A similar derivation for the time-dependent Hartree-Fock-Bogoliubov (TD- 
HFB) equations would be desirable for applications both in nuclear physics 
[2] and in the physics of trapped fermion droplets [3], however the Wigner- 
transformed TDHFB equations are a rather complicated set of coupled differen- 
tial equations [3] for the normal and anomalous densities p(r, p, t) and re(r, p, t) 
(p is real, while k is complex) and it is not obvious which approximations should 
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be introduced in order to obtain the equivalent of the Vlasov equation for su- 
perfluid systems. 

We have studied such an approximation that has the advantage of leading 
to simple solvable equations, even if it has the problem of violating particle- 
number conservation. We shall see that this problem can be solved within 
the model. Here we discuss this model for a one-dimensional system of spin- 
saturated fermions. Like for the normal Vlasov equation, understanding the 
simpler one-dimensional problem is an important preliminary step for study- 
ing three-dimensional spherical systems [4] . By comparing the response of the 
correlated system to an external perturbation with that of a normal system 
described by the ordinary Vlasov equation, we can learn about the effects of 
pairing correlations on the system dynamics. 

A fundamental role in determining the system eigenfrcqucncics is played 
by the boundary conditions imposed on the fluctuations of the density, here we 
make the comparison between correlated and uncorrelated systems by employing 
the simple fixed-surface boundary conditions of [4] , but fancier (moving-surface) 
boundary conditions might also be considered [5] (see also [6]). 



2 Static limit 

Bengtsson and Schuck [7] have studied the semiclassical limit of the static HFB 
equations (see also [1], p. 550). By following their method, we find the self- 
consistent equilibrium solutions 1 

Po(r,p) = -(l- E(r>p) ) (2) 
, x A (r,p) 

Ko(r ' p) = -2^ir (3) 

where 

• ho(r, p) is the self-consistent Hartrec-Fock equilibrium hamiltonian, 

• the chemical potential /j, is determined by the number of particles: 

9 



J drdpp (r,p) (4) 



(27T?i) 3 

(g is the number of fermions in a phase-space cell, g = 4 for nucleons), 
• E(r, p) is the quasiparticle energy, defined as 

E(r, p) - ^Ag(r,p) + (fco(r,p)-//) 2 , (5) 



x Note that our sign of kq differs from that given in Eq. (9) of [7], this means that our k 
corresponds to — k of Refs. [7], [3] and [1]. 
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• Ao(r, p) is the equilibrium pairing field, which is related to Ko(r,p) also 
by the self-consistency relation ([1], p. 550) 2 



A (r,p) 



1 



(2tt?i) 3 



dkv(\p - k|)K (r, k) . 



(6) 



We take a more phenomenological approach and replace the pairing field A (r, p) 
with the phenomenological pairing gap, in heavy nuclei A « lMeV. In this case 
the static solutions become 



Me) = i( 
«o(e) = 



1 



e - M 
E(e) ) ' 



me) ' 

£(e) = VA 2 + (e - A*) 2 , 
with e = /io( r 5 p) the particle energy. Note that 

E 2 (e)dp (e) 



«o(e) 



(7) 

(8) 
(9) 

(10) 



In Fig.l the functions po(c) and Ko(e) are plotted for typical values of pa- 
rameters (j, = 33MeV and A = lMeV. 
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Figure 1: Normal (solid, upper) and anomalous (solid, lower) densities as a 
function of particle energy, expressed in MeV. The dashed curve shows the 
normal distribution for A = 0.01 MeV. 



2 In view of the previous remark on the sign of k, here the sign of v(\p — k|) must be 
opposite to that used in [1] 
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3 Linear response 



We want to study the linear response of a system that, in equilibrium conditions, 
is described by the equations above. We assume that we have a finite, spin- 
saturated system of fcrmions at zero temperature with pairing correlations that, 
at time t = is subject to a weak perturbing external force generated by a 
driving field of the kind 5V ext (r, t) = (3S(t)Q(r) (the parameter (3 determines 
the intensity of the external force, the ^-function its time dependence and Q(r) 
its spatial distribution). We want to determine the fluctuations of the density 
of the system at point r as a function of time for t > 0. This will be done in 
linear approximation, that is, we neglect terms of second order or higher in the 
fluctuations. Thus we have: 

h(T,p,t) = h (r,p) + Sh(r,p,t)=e + SV ext (r,t), (11) 
p(r,p,t) = p {e)+Sp(r,p,t), (12) 
/c(r,p,i) = K (e) + SK r (r,p,t)+iSKi(r,p,t) (13) 

and, from the TDHFB equations, the following coupled equations for 5p(r, p, t) 
and diii(r, p, t) can be obtained: 

d t Sp(r,p,t) = {h,p}-2jSKi(r,p,t) (14) 



Sh(r,p,t)). (15) 



d t 6 Ki (r,p,t) = 2-^(-[<fp(r,p,i)+Jp(r,-p,i)] 

dp 
de 

These two coupled equations can be interpreted as an extension of the ordinary 
Vlasov equation of normal systems to systems with pairing. In the limit A — > 0, 
the first equation gives the usual Vlasov equation, while the second equation can 
be ignored. When A^O the two equations are coupled. Here we consider only 
the approximation Sh(r, p, t) = 5V ext (r, t), which corresponds to the single- 
particle approximation of the quantum approach, taking into account also the 
mean-field fluctuations induced by the external field leads to collective effects 
(see [4]). 

Taking the Fourier transform in time of Eqs. (14, 15), gives 

-iu6p(r,p,u) = {h,p} -2^-5Ki(r,p,uj) , (16) 



-iui6Ki(r, p,u;)) = 2-^-(-[5p(r,p,w)+5p(r,-p,u)] 
de 



^5h{v,p,oj)) , (17) 



or, for u! ^ 0, 
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-iw6p(r, p, w) + {5p, h } = -iwd 2 ^[Sp(r,p,u) +5p(r,-p,u)] 

- ^({ho,6h}-iLud 2 5h) , (18) 

with 

<"> 

and 

0(e) = 2^ . (20) 

In Eq. (18) we have used the relation {po,Sh} = ^-{h 0l Sh}. The equation 
corresponding to Eq.(18) for normal systems can be formally obtained from this 
equation by letting d 2 — > 0. (cf. Eq. (2.12) of [4]). Thus the effect of pairing 
in the present approximation is only that of adding the terms containing d 2 in 
Eq. (18). 



4 One— dimensional system 

For one-dimensional systems it is convenient to make the change of variables 

2 

(x,p x ) — ► {x, e), with e — ^ + V (x). Then Eq. (18) gives the following system 
of coupled differential equations: 



where 

/± - 6p(x, ±^2m[e-V (x)}, w) (23) 
and the inhomogeneous term 

B ± (a;) = B{x) ±d 2 C(x) (24) 

contains both the inhomogeneous term of the normal Vlasov equation 



B( 



(25) 



and an additional term proportional to the anomalous density (cf. Eq. (10)) 
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We have neglected the contribution of mean field fluctuations to the inhomo- 
geneous terms B(x) and C(x), so we limit our discussion to the zero-order 
approximation of [4]. In this case the coupled system (21, 22) can be solved 
exactly and its solution allows us to write an expression of the correlated prop- 
agator analogous to that obtained in [4] for the uncorrelated propagator. 
The uncorrelated propagator obtained in [4] was 

D°(x,x',u;) = ^ (27) 



dpo — 2nu> a cos[nuj T(x)] 1 cos[naj T(x')] 



de ^ T v(e,x) u — nu + is v(e,x') 

n 

while the correlated propagator given by the solution of the system (21, 22) is 

&>{x,a/,w) = JL (28) 

2 I d sr^ ~2&n cos{nLd T(x)] 1 cos[nuj T(x')] 
J de ^ T v(e,x) lu — u> n + is v(e,x') 

By comparing the expressions of the two propagators, we can immediately 
appreciate some of the changes induced by the pairing correlations: 

• The equilibrium distribution is changed from the familiar Fermi-gas step 
function (dashed curve in Fig. 1) to the smoother curve also shown in 
Fig. 1, so the function which is a <5-function in the uncorrelated case, 
becomes also smoother and this smears the response. 

• The position of the propagator's poles are determined by the eigenfrequen- 
cies of the system. For the uncorrelated propagator they are u) n — ncj , 
where n is an arbitrary integer and o>o = 2ir/T is determined by the system 
size (for given e), while for the correlated propagator, the cigenfrcqucncies 
become 

in = (~~) 2 for n ^0, (29) 

u>„ = ±fi(e) forn = 0. (30) 

A closer analysis of the correlated response function 

S°(huj) = ~lm j dx j dx'Q(x)D°(x,x',hLj)Q(x') (31) 

reveals other interesting properties of the correlated system: 

• the correlated response function displays a gap of ~ 2A at low excitation 
energy; 

• the lack of self-consistency due to the constant-A approximation gener- 
ates some spurious strength in the correlated response. 
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There are two sources of spurious strength: one is related to the mode n = 
and its possible contribution to particle-number non conservation. However, 
even if the contribution of this mode is omitted from the propagator (28), the 
resulting strength function gives an energy-weighted sum rule 

r°° a 
Mi = / dhLuhLuS°(huj) w Mi[l + 14.66(- — ) 2 ] (32) 
Jo 

which differs from the uncorrelated sum rule 

9 AT 2 

Ml = U?^. (33) 
3 m 

What is even more disturbing is the fact that the enhancement factor depends on 
the size of the system. Both problems can be solved by introducing a modified 
propagator 

Dl hys {x,x',uj) = ^ (34) 

dpo y^.ntJo 2 — cos[nuJoT(x)] 1 cos[na>o T (a ; ')] 

de ^ lu„ T v(e,x) io—u n + ie v(e,x') 

This modified propagator contains no contribution from the mode n = and 
this ensures particle-number conservation, since 

5A(lo) = j dxdp x 5p(x,px,u)) = J dxdx' D® hys (x,x' ,u)Q(x') (35) 

vanishes for any external field Q(x') if only modes with n / are included in the 
propagator. Moreover, it can be easily checked that the strength function given 
by this propagator satisfies the same sum rule as that given by the uncorrelated 
propagator. 

The following two figures show examples of correlated and uncorrelated re- 
sponse function for two systems of different size. The comparison is made by 
assuming hat the systems have the same value of fi, hence the number of parti- 
cles is slightly different. 

Figure 2 shows a comparison of the first peak in the correlated and uncor- 
related strength functions for the external field Q(x) — x 2 and an equilibrium 
mean field Vo(x) of the square-well type. The response functions have been 
divided by a factor -^j-L 4 , where L is the size of the square-well mean field, 
which has been chosen so that TiwoGu) = lOMeV. 

Figure 3 shows another example of response to the same external field Q(x) = 
x 2 , but in this case the size of the potential well has been chosen so that Tiuj {p) = 
1 MeV. All other parameters are unchanged. In this case ?«i„ = i(/i) w 2.23 MeV 
and the main peak appearing in this figure (solid curve) corresponds to the 
mode n = 1, the smaller bumps at larger energy are due to the modes with 
larger values of n. Note the gap of ~ 2A in the low-energy region, the small 
tail extending into the gap is due to the finite value of e used in the calculation. 
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Figure 2: Uncorrelated and correlated response functions (renormalized), as 
a function of excitation energy Tilo expressed in MeV. The peaks shown here 
correspond to the mode n = 1 of the two propagators (27) and (28). The peak 
on the left should be a <5-function and it represents the response of a normal 
system given by the propagator (27); for numerical reasons we have used a finite 
value of e = 0.1 MeV, thus inducing an artificial smearing of the response. The 
peak on the right is obtained from the propagator (34) with the same value of 
e. Its larger width is an effect of the pairing correlations. There are other similar 
peaks around Tilo — 20, 30 . . . MeV, but their strength decreases rapidly with 
increasing n. 
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Figure 3: Same as Fig. 2, but with Hu)o((j,) = 1 McV and including modes up 
to n = 10. The main peak of the uncorrelated strength function (dashed) at 
hw = 1 McV is pushed to higher energy by the pairing correlations and a gap 
of about 2A is created at low energy. 
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